Thermophoresis of Brownian particles driven by coloured noise 
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Brownian motion of microscopic particles is driven by collisions with surrounding fluid molecules. 
The resulting noise is not white, but coloured, due, e.g., to the presence of hydrodynamic mem- 
ory. The noise characteristic time-scale is typically of the same order of magnitude as the inertial 
time-scale over which the particle's kinetic energy is lost due to friction. We demonstrate theoreti- 
cally that, in the presence of a temperature gradient, the interplay between these two characteristic 
time-scales can have measurable consequences on the particle's long-time behaviour. Using ho- 
mogenization theory, we analyse the infinitesimal generator of the stochastic differential equation 
describing the system in the limit where the two time-scales are taken to zero keeping their ra- 
tio constant and derive the thermophoretic transport coefficient, which, we find, can vary in both 
magnitude and sign, as observed in experiments. Studying the long-term stationary particle dis- 
tribution, we show that particles accumulate towards the colder (positive thermophoresis) or the 
hotter (negative thermophoresis) regions depending on their physical parameters. 

I. INTRODUCTION 

A microscopic or nanoscopic object immersed in a fluid, e.g., a Brownian particle or a biomolcculc, undergoes a 
permanent thermal motion. This motion is the result of the collisions with the fluid's molecules and is typically 
modelled as driven by a white Gaussian noise However, this driving noise is actually coloured, i.e. it has a 
characteristic non-zero correlation time r, on a very short time-scale of the order of tens of nanoseconds, due, e.g., 
to the presence of hydrodynamic memory This time-scale is similar to the particle's inertial relaxation time, i.e. 
the characteristic time for loss of kinetic energy through friction 17 = 771/7, where m is the mass of the particle and 7 
the friction coefficient Q. As we will theoretically demonstrate, the interplay between these effects, despite occurring 
on time-scales that might not be in themselves directly accessible experimentally, can have measurable effects on the 
particle's long-time behaviour. 

In this letter, we consider the dynamics of a Brownian particle driven by a coloured noise when it is immersed 
in a fluid where a temperature gradient is present. We find that due to the interplay between its two characteristic 
time-scales, r and cr, the particle can exhibit a directed motion in response to the temperature gradient; furthermore, 
studying the long-term stationary particle distribution, we find that particles can accumulate towards the colder 
(positive thermophoresis) or the hotter (negative thermophoresis) regions depending on their physical parameters 
and, in particular, on the dependence of their mobility on the temperature. We remark that, as demonstrated by the 
mathematical analysis presented below, the presence of a coloured noise, as opposed to a white noise, is crucial for 
the emergence of such thermophoretic effects. The velocity of this motion can vary both in magnitude and sign, as 
observed in experiments 0-0]. As a result, the particle moves either towards the cold or the hot side, quite similarly 
to what happens in the presence of an external driving force, e.g., gravity or electric fields; however, in this case, no 
external force is actually acting on the particles Q. 

II. MATHEMATICAL MODEL 

For a spherical particle of radius R immersed in a fluid of viscosity /i, which in general depends on the absolute 
temperature T, i.e. \i = fJ-{T), the friction coefficient 7 satisfies Stokes law, 

7 (T) = 6nn(T)R, (I) 

and the diffusion coefficient D is related to 7 by the fluctuation-dissipation relation 0, 

D(T) = (2) 

For ease of argument, we will assume that the particle's motion is one-dimensional in a horizontal direction perpen- 
dicular to gravity, with position denoted by Xt € R for all times t > 0. We also assume viscosity to depend only 
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on temperature, and fluid thermal expansion and convection to be negligible. If more complex models are needed, 
e.g. to account for interactions between the fluid and the particle or for the thermal expansion of the fluid, then 
the above model may need to be modified and the results may change, but the approach to address such problem 
will be the same as described in this Letter. The relation in eq. ^ assumes local thermodynamic equilibrium [9|, 
which implicates that the temperature gradients should not be too steep; such conditions have been shown to be 
experimentally verified, e.g., in Ref. The resulting motion of the particle is governed by the stochastic Newton 
equation: 

Vt 



mx t = -y(xt)x t + j(x t )\/2D(xt)—j=, (3) 



with initial conditions xo — x* and xq — v*. The coloured noise r\t is an Ornstein-Uhlenbeck process (OUP) defined 
by the stochastic differential equation (SDE): 

dr)t = --Vtdt + J-dW t . (4) 

T V T 

where r > is the noise correlation time and Wt is a standard Wiener process. Its stationary solution is a zero-mean 
Gaussian process with E[(rj t r] s )/T] = (1/r) exp { — -\t — s|} and its covariance function converges to the delta function 
as r tends to zero. 

Typically, the hydrodynamic and inertial memory time-scales, i.e r and a, are very fast - in particular, much faster 
than the typical time resolution at which the particle Brownian motion is experimentally sampled - and their effects 
accumulate over the longer diffusive time-scale and eventually results in transport dynamics. Therefore, following 
Langcvin's approach, it is customary to drop the inertial term, i.e. set the left hand side to zero in eq. <|3j> PJ and take 
the time correlation of the noise to zero. However, the limit of eq. ([3]) as m — > 0, i.e. a — > 0, has to be studied with 
care, requiring a nontrivial computation, as, in general, similar limits involve additional drift terms pdl4l3| . Here we 
are interested in the long-time behaviour of x t and in how the particle undergoes a deterministic drift in response 
to a temperature gradient as both the inertial and the noise characteristic times are taken to zero. To this end, we 
will study eq. ([3]) analysing the convergence of the infinitesimal operators of the corresponding diffusion processes 
(backward Kolmogorov equations), thus following the well-known methods from homogenization theory As 
a — > and r — > 0, we expect the limiting equation to be of the form: 

dx t = A(x t ) dt + B(x t ) dW t , (5) 

where the effective coefficients A(x) and B(x) depend in general, as shown below, on the relative rate at which a and 
t go to zero. Therefore, we define 

* = - = (6) 
r 7T 

which is a dimensionless parameter and, in general, depends on x, i.e. 9 = 9(x). We note that, for consistency, 
we choose to write the limiting eq. ([5]) with the ltd convention - a different choice would change the coefficients in 
the following discussion, but not the physical meaning. The effective drift term A(x) reflects contributions resulting 
from the coloured nature of the noise, as well as from the small mass approximation. We stress that the particle's 
motion is perpendicular to gravity and A(x) represents a real drift in the particle position, which arises despite the 
absence of external forces acting on it. This is a consequence of the system not being in thermodynamic equilibrium, 
as manifested by the presence of a temperature gradient 0]. 

For constant viscosity, i.e. fJ-(x) = [Iq and thus "f(x) = 70, the simultaneous limit as r — > and a — > was studied 
by Freidlin [f2j|: by taking t — > first and then a — > 0, he obtained the zero-mass limit eq. (0 with A(x) = 0, 
while, by taking a — > first and then r — > 0, A(x) = B ^ x '^ ( x ) . Subsequently, Kupferman, Pavliotis and Stuart (Tl| 
showed that, if in the limit 9 is constant in x, all the intermediate cases can be obtained as 9 is varied. In a previous 
article [HI], we studied the limiting equation for a general class of equations similar to eq. ([3]) but driven by a white 
noise. In this letter, we study the analogous problem for the coloured noise case. 
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III. MULTISCALE ANALYSIS 



To simplify the following analysis, we rewrite eq. ([3]) as a set of first order SDEs, setting m = ■yOr [eq. ([6])] and 
substituting u[ - 



dx\ 



duT 



drjt 



uT dt 



-$m dt 



2D(xl) 1 (xp ru 
9 T 



dt 



(7) 



4 dW t 



with the initial conditions 



v* and 770 a normal random variable with mean zero and variance - 
independent of the Wiener process Wt ■ The notation x[ is introduced to distinguish the full 3-dimensional SDE with 
coloured noise and mass, from xt which, from now on, represents the limit as t, <j — > 0. 

To determine the coefficients A(x) and B(x) for the limiting eq. we use a multiscalc analysis of the backward 
Kolmogorov equation associated with eqs. ((?]) fl4T - [l6| : we will first write the backward Kolmogorov equation for 
cqs. ([7]); wc will then derive its limit as r — > keeping constant; and, finally, we will determine A(x) and B(x) from 
this limiting Kolmogorov equation. 

Let g(x', it', 77', t'\x, u, 77, t) be the probability density of the distribution of the position (x'), the rescaled velocity 
(u') and the coloured noise (77') of the particle at time t' given their values (x, 14,77) at a time t < t' . The backward 
Kolmogorov equation for eqs. ([7]) is: 



dg 
dt 



—^o H — i=Li 

t Jr 



(8) 



with Lr, = 



2D j „ 



^■ — 277^+2^7 and ii = -^=^. We will take advantage of the fact that eq. ([5} involves 

derivatives with respect to the x,u,r},t variables only and shorten the notation by writing g(x', u', 77', t'\x, u, f], t) = 
g(x,u,rj,t). Furthermore, we will not always explicitly indicate the dependence of 7, D and 8 on x. 

Following the general multiscale analysis ansatz |14l - fl6l ] , we postulate the solution to the Kolmogorov equation as 
a sum of an asymptotic series, 



g = go + V^gi + rg 2 + 

and we substitute it in eq. ([5]). By matching corresponding powers of \fr, we obtain the following equations: 

^oSo = 0, 
Lo<?i = —Ligo, 

dgo T . T 
at 

Eq. (|10p implies that go(x,u,r]) = go(x), since all derivatives of Lq act on u and rj only. 
A solution of eq. (fTTj) has the form, 



where $ solves the cell problem (see [l6| sec. 11.3) defined as 

- Lo$(x,«, 77) : 

We look for a solution of the form 



dgo 

dx ' 



70 



-u + V2Dh(rj) 

1 



(9) 

(10) 
(11) 
(12) 



(13) 



(14) 



(15) 



Substituting eq. (fT5|) into eq. (p^4|) . we get the equation + 2h vv = 77 with a solution h = 77/2 and, therefore, 
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$(X,U,7]) = J-U 



-'/■ 



(16) 
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The solvability condition for cq. (|12l) follows from the Fredholm alternative [16H18I ]: 



(17) 



for all p such that L^p = 0, where Lq is the dual operator of L$. Thus, we solve the dual problem L^p = with the 
dual operator dchncd as 



t " 



d_ 

du 



2Dj 



' drj y 



This problem is equivalent to finding the stationary distribution for a two dimensional OUP in u and 77, with 1 as a 
parameter, i.e. a zero-mean Gaussian probability density p(u,rj;x) = C(x) exp { — \si±u 2 — s^urj — ^S22?7 2 }, where 

1 ( l+29 \ 2 _ 1 1+29 J 

Wj \-28~) ' Sl2 - and 



C(x) is a normalizing factor. Solving for the parameters s+j we obtain: sn 



V2D 7 6/ 26 

s 22 = ^jf-- Using this expression for p we can evaluate the integral in eq. (fT7|> and obtain the backward Kolmogorov 
equation for 170 



dt - A{X) dx +B {X) dx 2 ■ 



(18) 



where, A(x) = J* R2 ^^^p(u,n; x) dud-q and B{x) 2 = 2 j" R2 -^=$ / o(it, 77; x) dudij. Finally, we can compute the 

covariances with respect to the stationary density p, i.e. E[urj\ = ^^^e ail< ^ E[u 2 ] = and obtain explicit 

formulae for the coefficients of the effective SDE ([5]): 



A(x) 



D'(x)j(x) - 46>(x)£)(x)7 / (x) 
2 7 (x)(l + 26»(x)) 



(19) 



and 



B(x) 2 = 2D{x). 



(20) 



We remark that the diffusion term in the limiting eq. is the same as in the initial eq. ([3]). The precise meaning 
of the limiting eq. (|5|) is that the solution x\ of eq. (|7|) converges in law to the solution x t of eq. ([5]) with the same 
initial condition xo = x* as r, a — > 0. A rigorous proof is justified by techniques in [l9j . 



IV. PHYSICAL INTERPRETATION: DRIFT AND PROBABILITY DENSITY 

Since the effective limiting SDE [eq.©] has been constructed using a stochastic integral with the Ito convention, 
the expected position is 



E[x t ] = E[x ] + E 



A(x s ) ds 



(21) 



The sign of A(x) determines the direction towards which a Brownian particle is expected to travel. Therefore, if, e.g., 
we have A(x) > 0, the particle will on average travel towards increasing x until it reaches some boundaries, which 
can be either absorbing boundaries or reflecting boundaries. 

In the presence of absorbing boundaries, particles disappear from the system as soon as they reach a boundary. 
Therefore, the sign of A(x) determines the boundary at which particles arc preferentially absorbed, as schematically 
shown in figure [TJa). This situation can be experimentally realized, for example, within a relatively long microfluidic 
channel where particles are steadily injected at a certain position, e.g., x = 0, and removed once they reach either end 
of the channel. In the presence of a temperature gradient, these particles move towards increasing or decreasing x 
depending on the sign of A(x) and, therefore, can be sorted and classified on the basis of their physical and chemical 
properties that influence A(x). Such a sorting process is more efficient the longer the channel and the smaller the 
noise term, i.e. B(x). 

In the presence of reflecting boundaries, particles are reflected back into the system when they reach a boundary. 
In this way, once a particle has interacted with the boundaries multiple times, the particle's position xt reaches a 
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FIG. 1: (a) Schematic representation of the behaviour of thermophoretic particles in the presence of absorbing boundaries: 
black (white) particles have Ai,(x) > (A w (x) < 0) [eq. (|19|l ] and are mostly pushed towards the right (left) boundary, where 
they are removed from the channel, (b) Schematic representation of the behaviour of thermophoretic particles in the presence 
of reflecting boundaries: the particles eventually reach a steady state probability density poo{x) [eq.l|24|)]. which is different for 
particles with Ab(x) > (black) and A w (x) < (white). 



steady state probability density p oa (x), as schematically indicated in figure QJb). The time-dependent p(x,t) is the 
solution to the following Fokkcr-Planck equation (also known as forward Kolmogorov equation) , 



dp 
dt 



d d 2 
- = --[A(x)p] + w 



B{xf 



-P 



(22) 



with a given initial condition p(x, 0) = po(x). Since the Fokker-Planck equation is deterministic, its solution, i.e. the 
evolution of the probability density over time, does not involve any randomness. As t — > +oo, the solution to eq. (|22j) 
converges to the steady state distribution p oc (x), under certain conditions [2(| sec. 5.2]. If we assume the motion of 
the particle to be restricted to the interval (a, b), a < b, then we can solve the stationary Fokker-Planck equation for 
the SDE ([5]) [2l[. Given A(x) and B(x) the stationary solution is 



Poo 0*0 



C 



B( x y 



■ exp 



B( X y 



dx 



which can also be expressed as a function of D(x), j(x) and 9(x) as 



p 00 (x) = CD(x) 2 + 4 »<->7(a;) 



(23) 



(24) 



where C is a normalizing constant. This situation can be experimentally realized in closed systems where a temperature 
gradient is present. Interestingly, this is the case of most experiments performed to study thcrmophoresis and the 
Soret effect, where a suspension of particles in a thermal gradient is given enough time to relax to its steady state 
distribution @, 0] . 



V. RESULTS AND DISCUSSION 



For the following discussion, we will express the thermophoretic drift as a function of T and p(T), and since p(T) 
is interpreted as an expansion around some temperature To, we will write AT = T(x) — To, and -^p(T) = //(T)AT. 
From eq. (|TT))) . using cqs. (Q]) and we find that the effective thermophoretic drift is 

_ p(T)-p'(T)AT[l + A9] 
A(x)-k B l 12nR ^ {T)[1 + 20] > I 25 ) 

where we have set T = T(x) and 9 — 9(x), suppressing the dependence on x for brevity and p'(T) and T' are 
derivatives with respect to T and x respectively. Eq. (|25[) shows that the thermophoretic drift is determined not only 
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FIG. 2: /i(T) constant, (a) Thermophoretic drift A(x) for various 6 (dashed, dashed-dotted and dotted lines) in the presence of 
a temperature gradient (solid line), (b) Corresponding theoretical steady state distribution (lines) are always peaked towards 
the cold side. The symbols represent the steady states distributions resulting from Brownian dynamics simulations. 



by T', but also by the dependence of /z on T. Interestingly, if fi'(T) > and /Lt(T) > n'(T)AT, from eq. ([25]). there is 
a critical 8 denoted 8 C such that the drift A(x) changes sign: 



»(T)-»'(T)AT 
V(T)AT 



(26) 



The stationary density is given as 



Poo (x) = C 



T(a;) 1 + 4 « 



(27) 



which is an inverse power of T unless /Lt(T) is at least quadratic in T. 

Below, we consider in more detail the three simplest cases, i.e. /i constant, /i linear in T and /.i quadratic in T. 
We stress that, even though these are the first three orders of approximation, it might be necessary to consider more 
complex situations in real applications because of the large range over which T can vary. 



A. n(T) constant 

In the simplest case /j, does not depend on T, i.e. jit(T) = /Ltn. > 0, and the thermophoretic drift term (|25|) becomes 

A(x) = k -^- -. (28) 

v ; 127rii/xo(l + 26) v ' 

We stress that in this case 8 is independent of x. For example, let us consider the temperature gradient [grey solid 
line in figure |2{a)]. The associated thermophoretic drift [ea.(|28p] is presented in figure a) for three values of 8, 
corresponding to the situations where the time correlation of the coloured noise dominates (8 = 0, blue dashed 
line), the time-scale of the particle inertia dominates (8 = +oo, red dotted dashed line) and the two time-scales are 
comparable (0 = 1, green dotted line). For all 8 > the drift term has the same sign as T", thus imposing a drift 
on the Brownian particle instantaneous flow towards the hotter region. Equation (|25[) features two terms, one term 
dependent on the frictional gradient and the second on the coloured noise; for \i constant, since there is no frictional 
gradient, the drift is only driven by coloured noise and the resulting flow increases with decreasing 8. 

If now the particles are allowed to interact repeatedly with the boundaries, they will eventually reach their stationary 
distribution, using eq. ([2^)l , 

1+49 

Poo (x) (xT(x)-^. (29) 

As shown in figure [5Jb), the particles will accumulate towards the areas of low temperature for all 8 > 0, which 
is in agreement with most experiments 0. This result is in striking contrast with the fact that the instantaneous 
thermophoretic drift actually pushes the particles in the opposite direction towards the hotter regions. Such a 
difference between the instantaneous drift and the long-term stationary distribution has also been observed in systems 
at thermodynamic equilibrium [22 25|. 
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FIG. 3: Same as figure [2] for fi(T) linear. Note in (a) the sign change of A(x) as 6 crosses C 
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FIG. 4: Same as figure [2] for n(T) quadratic. Note in (b) the change of the distribution peak from cold to hot as a function of 



B. fj,(T) linear 

We now consider the case when //(T) is linear, i.e. fj,(T(x)) = fj,o + ^\AT(x) > 0. Interestingly, if /^i > (e.g. Ref. 
Hfl), A(x) changes sign at a critical value of 9 [eq. (p6)) ] 



Mo 



(30) 



The resulting effective thcrmophoretic drift (p25)) is shown in figure El(a) for the cases of 9 = < 9 C , 9 = 1 > 6 C , and 
9 = co. In particular, 



and 



A(x) = 



A(x) = 



127Ti?(^ + MiAT) 2 



for 9 = 



67ri?(^ + MiAT) 5 



for 9 



-co, 



(31) 



(32) 



whose dependence on the thermal gradient clearly shows opposite signs if p,\ > 0. 
For the stationary distribution [figure GDJb)], 



Ax>(x)cx \fi AT(x)- {1+4e) +fi 1 AT(x) 



-46 



(33) 



For \i\ > it is clear that will be an inverse power of T(x). If fii < 0, since 7(2;) > for all x, then (Iq > \fiiAT(x)\ 
and /j,qAT(x) ( 1+4e ^ > \fx 1 AT(x) Thus poo will be an inverse power of temperature for all admissible /io,A*i- 
This suggests that the particle will more likely be found in the colder regions. Interestingly, unlike the expected drift, 
derived from eqs. (|31[) - (|32[) , there is no qualitative change in behaviour at 9 C . 
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C. fJ>(T) quadratic 



We finally consider the case when fi(T) is quadratic in T, i.e. /i(T(x)) = Mo + MiAT(x) + M2AT(x) 2 . Again, A(x) 
changes sign at [cq. (126)) ] 



Mo - M2AT 



(34) 



4(/i!AT + 2^Ar) 

if fi > /i 2 AT 2 and MiAT + 2yU 2 AT 2 > 0, or Mo < M2AT 2 and Mi AT + 2^ 2 AT 2 < 0. The resulting effective 
thermophoretic drift (j2"5)l for the extreme cases is 

= fcgT ( Mo - m 2 AT 2 ) for = 

127ri?(/i + /iiAT + M2 AT 2 ) 2 

and 

am k B T'^ 1 AT + 2^ 2 AT 2 ) 

A(x) = 5 — tor ft — > +00, (36) 

V ; 67ri?(Mo + mAT + M 2 AT 2 ) 2 V ' 

whose dependence on the thermal gradient shows opposite signs for mo > M2AT 2 . In figure EJa), we study the case 
Mo < M2AT 2 and mi AT + 2m 2 AT 2 . 



The stationary distribution is 



Poo oc 



MoAT -(i+40) + MlAT -4e + M2AT (i 



(37) 



where it is clear that there are combinations of mo> Mi an d M2 for which the density will incur a transition from 
peaking at colder regions to hotter, as shown in figure BJb). This may lead to interesting behaviours as a function of 
the temperature. For example, in reference [27|, DNA particles (8 <C 1) change from accumulating in colder regions 
when the minimum temperature is T = 276 K to accumulating in warmer regions when a new minimum temperature 
T = 293 K and the gradient is left unchanged (2 K between the colder and warmer regions). In agreement with this 
experiment we predict, for T = 276 K, since, expanding fi(T) around T = 273, Mo + Mi AT > M2AT 2 , the stationary 
distribution Moo [cq. (|37[) ] peaks in colder regions for all 8, and, for T = 293 K, since M2AT 2 > Mo + MiAT, the 
stationary distribution peaks in hotter regions for 8 < 1/4, which is verified in this case. However, we note that the 
model in this paper leads to a different explanation than the one in reference [27l ] and there may be numerous other 
factors influencing thcrmophoresis including thermal expansion of the fluid and convection. 



VI. CONCLUSIONS AND OUTLOOK 



In this letter, we have given a systematic analysis of a system with two short time-scales: the time correlation of 
the coloured noise t and the inertial relaxation time a. We have derived the effective thermophoretic drift A(x) and 
the steady-state probability distribution Moo (2) in the limit as these time-scales go to zero. We have finally applied 
these results to study the thermophoretic motion of a Brownian particle in a temperature gradient, showing how, in 
agreement with experiments, Moo (2) tends to be peaked towards the colder region (positive thcrmophoresis), but can 
switch to hotter regions (negative thcrmophoresis) under the right conditions. As possible future lines of research, 
also noises that are not OUP can be considered: these might be particularly promising to model the effects of the 
chemical interaction between a particle and a solvent, of the viscoelasticity of the medium or of depletion forces due, 
e.g., to polymers or small particles in the solution. 
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